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We point out that the non-adiabatic orientation of quantum rotors, produced by ultrashort laser 
pulses, is remarkably enhanced by introducing dipolar interaction between the rotors. This enhanced 
orientation of quantum rotors is in contrast with the behavior of classical paired rotors, in which 
dipolar interactions prevent the orientation of the rotors. We demonstrate also that a specially 
designed sequence of pulses can most efficiently enhances the orientation of quantum paired rotors. 
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I. INTRODUCTION 



Considerable attention has been paid to the ability of 
intense laser fields to align or orient polar molecules; in 
such fields, the molecules experience torque arising from 
the dipolar interaction with electric fields [1]. One ap- 
proach to achieve adiabatic molecular alignment is to use 
a nanosecond-pulse laser [2-4]. While adiabatic align- 
ment disappears after the pulse is turned off, ultrashort 
pulses (100 fs or less) can excite rotational wave pack- 
ets of quantum rotors, thus yielding a noticeably aligned 
shape after the pulse is off [5-11]. Furthermore, a spe- 
cially designed sequence of pulses is known to achieve an 
enhanced angular focusing in quantum rotors [12, 13]; 
this has been realized experimentally in optical lattices 
[14]. These types of alignment, which differ from the adi- 
abatic, are important for manifold applications requiring 
transient molecular alignment under field-free conditions, 
such as the generation of laser pulses [15, 16] and the con- 
trol of high harmonic generation as a source of coherent 
radiation [17]. 

In most studies of the molecular alignment, dipolar 
interaction between polar molecules has been neglected. 
In this case, the rotational dynamics of polar molecules 
can be simply analyzed by using an isolated kicked-rotor 
model [18]. The quantum kicked rotor and its classical 
analog have long served as a paradigm for quantum and 
classical chaos. In contrast, the quantum dynamics of in- 
teracting kicked rotors has been less well studied so far, 
though intriguing phenomena of interacting quantum ro- 
tors are expected to emerge. Very recently, a study of the 
center-of-mass motion in two coupled kicked rotors has 
revealed that the decoherence effect induced by the in- 
ternal degree of freedom enhances the quantum-classical 
correspondence in the dynamics of rotors [19]. In addi- 
tion, anomalous dielectric responses have been pointed 
out in two coupled dipolar rotors [20] . These results im- 
ply that coupled kicked rotors, even in only two rotors, 
exhibit peculiar behaviors different from those of isolated 
kicked rotors. 

In the present work, we theoretically investigate the 
quantum dynamics of coupled (paired) kicked rotors sub- 
ject to an ultrashort laser pulse [5- function kick). We 
find that the dipolar interaction between rotors remark- 



ably enhances the transient orientation of paired rotors 
produced by a <5-kick. This enhancement of the orienta- 
tion of quantum rotors contrasts with the results with the 
classical treatment for paired rotors; in the latter case, 
the dipolar interaction inevitably hinders the orientation 
of coupled rotors. Furthermore, we demonstrate that the 
orientation of paired rotors can be further enhanced by 
applying the accumulative squeezing scheme proposed in 
Rcf. [12]. Our findings enlighten the study of transient 
orientation of interacting polar molecules. 

This paper is organized as follows. Section II describes 
the Hamiltonian of the paired-rotor system together with 
its analytical solutions for eigenenergies and their eigen- 
functions. The time development of paired rotors for 
post-kicked times is given in this section. Section III an- 
alyzes the time dependence of the orientation factor in 
both quantum and classical paired rotors. The physical 
origin of the enhanced orientation in quantum rotors is 
discussed on the basis of the spatial profile of the prob- 
ability density of paired-rotor wavefunctions. Section IV 
gives the conclusion. The paper contains two Appendices 
with the details of the calculations. 



II. PAIRED-ROTOR SYSTEMS 

A. The Hamiltonian 

Suppose that two rotors carrying dipole moments \x are 
arranged as shown in Fig. 1. For arrangement (a), both 
rotors rotate in a plane, while in the case of (b), both ro- 
tors belong to an identical rotation axis [21]. The Hamil- 
tonian for the system is given by H = H i + H 2 + W\2 , 
where Hi is the Hamiltonian for the i-th kicked rotor and 
W±2 represents dipolar interaction between rotors. The 
term Hi can be written as 



Hi 
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V(fii,t), 



(1) 



where Li is the angular momentum operator and / is the 
moment of inertia of the rotor. When rotors are driven 
by a linearly polarized field, we can set 



(2) 



2 




FIG. 1: Definitions of arrangements of two dipolar rotors: 
(a) Two dipoles rotate in an identical plane, (b) Two rotors 
belong to an identical rotation axis. In both cases, the rotors 
are separated by R and interact via the dipolar interaction 
Wi2- The direction of the electric field E(t) is defined as 
shown. 



where E(t) is the field amplitude of short laser pulses. 
The direction of the field is fixed as shown in Fig. 1, 
causing the angular focusing of paired rotors at Ox — 
62 = 0. Assuming that the rotational radius of rotors is 
sufficiently smaller than the separation R between rotors, 
the dipolar interaction W12 is expressed by 
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(3) 



with the dielectric constant s. The vector R connects two 
rotational centers as denoted in Fig. 1. For simplicity, we 
rewrite the term (3) as 



W l2 

with the definition, 



Ed • F{0i,0 2 ) 



Ed — 



471-ei? 3 ' 



(4) 



(5) 



The quantity Ed determines the magnitude of the dipo- 
lar interaction between rotors, and plays a key role in 
the dynamics of coupled rotors, as we discuss later. The 
explicit form of F (6 1,62) is obtained straightforwardly 
from Fig. 1. The arrangement (a) gives the form, 

F(0i,6 2 ) = cos6»iCosfl 2 - 2 sin 9i sin 9 2 , (6) 

while, for arrangement (b), we have 

F '{Oit 02) — cos ^1 cos 2 + sin ^1 sin ^2 

= 008(0! - 2 ). (7) 



B. Eigenenergies and their eigenfunctions 

In the absence of the field E(t), eigenstates of paired 
rotors are analytically obtained by transforming variables 
into 2 )I2 and r\={Q x - 2 )/2. Substituting 

them into Eqs. (1) and (4), the Hamiltonian H is sepa- 
rated as [20], 



H = H~£ + H v , 



(8) 



Ek&_ 
2 do? 



E D c a cos2(a + a ); a = £,rj. (9) 



Here the quantity Ek = h 2 /(2I) represents the ki- 
netic energy. The parameters (c^, c v , £ , Vo ) equal 
(3/2,l/2,0,7r/2) for arrangement (a) and (0,1,0,0) for 
(b) . The separability of the Hamiltonian H allows us to 
write the paired-rotor wavefunction in the form rj) = 
tptiOPriiv)- Consequently the Schrodinger equation of 
paired rotors -/?<?■(£, 77) = E^{^,rj) can be decomposed 
into two independent eigenvalue equations expressed by 



da 2 



[e a - 2w a cos2(a + a )} if a = 0; a = £,77, (10) 



where e a — 2E a /Ex and v a = c a ED / Ek- The solu- 
tion of Eq.(lO) is given by the Mathieu function [22], 
whose explicit forms are given in Appendix A. The 
eigenenergies E of paired rotors are thus expressed by 
E=(e x + e y )E K /2. 



C. Time development of wavefunctions 

Let us consider the time development of the wave- 
function in paired-rotor systems after a <5-function kick 
at t = t. The wavefunction W(0i,02,t + ) immediately 
after the kick is related to that just before the kick, 
&(0i,02,t~), with a phase determined by 



<^(01,02,T + ) 



exp 



h 



{V(0 1 ,t) + V(0 2 ,t)}dt 



Substituting the definition (2) into Eq. (11) and trans- 
forming variables (#1, #2) — ► (£, rj), we obtain 

V(Z,V,t + )= eJn ( 1T C0S V exp(-mr ? )^,r ? ,r-), 

(12) 

where J n (z) is the Bessel function of n-th order. The 
quantity, 



P = 



fj,E(t)dt, 



(13) 



represents the strength of the pulse. In actual calcula- 
tions for Eq. (12), the summation of n can be truncated 
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at the finite value ±n c , since the magnitude of the Bessel 
function J n {^S- cos£) rapidly decays with increasing \n\. 
The time development of paired rotors for posi-kicked 
times is described by 



+i) = cxp 



:{H t + H v )t 



m,V,r + ). (14) 



The right-hand side of Eq. (14) can be analytically calcu- 
lated by expanding the function \P(t;, 77, r + ) by the Math- 
ieu function. The details of the calculation are presented 
in Appendix B. In the following, we set t = 0, and take 
Ek and K/Ek as units of energy and time, respectively. 
The initial state just before the kick is fixed in the ground 
state. 



III. ORIENTATION OF PAIRED ROTORS 

A. The orientation factor 

The degree of orientation for paired rotors is character- 
ized by the orientation factor 0(t) = (2 — cos 6\ — cos 62), 
where the angular bracket indicates to take the expecta- 
tion value. The factor 0(t) tends to be zero when the ori- 
entation of the rotors in the field direction of E is perfect. 
On the other hand, the factor equals 2 when the ampli- 
tude of wavcfunction 1 1^(6*1 , ^2 ) | is uniformly distributed 
in the 61-62 space. The strength of the dipolar interaction 
is characterized by the parameter T = Ed/Ek 1 i.e., the 
ratio of the interaction energy Ed to the kinetic energy 
Ek- 

Figure 2 (a) and (b) show the time dependence for the 
orientation factor 0(t) for paired rotors fixing the kick 
strength P — 10. Indices (a) and (b) in Fig. 2 correspond 
to arrangement of rotors (a) and (b) in Fig. 1. In both 
(a) and (b), the left figure plots the orientation factor 
for a long time scale, < t < 7.0, while the right one 
does so for a short time scale, < t < 0.15. We first 
discuss the two figures on the left, where two values of T 
are taken; the dotted lines display the orientation factor 
for r = 0, and the solid lines display that for T = 30.0. 
For r = 0, the time dependence for 0{t) for arrangement 
(a) is identical to that for (b), because the two rotors are 
no longer correlated via dipolar interaction. In this case, 
the orientation factor yields a simple form [23] of 



0(t) = 2 - 2Ji (2Psin£) 



(15) 



with the lowest value of 0(t c ) = 0.836 at the focal time 
t = f c = 9.2x 10~ 2 . 

For finite T's, the orientation factor exhibits somewhat 
complicated behavior different from that for isolated ro- 
tors. The left two figures in Fig. 2 exhibit the differ- 
ence of the time dependence for 0{t) between the case 
of T = 30.0 and that of Y — 0. For arrangement (a), the 
magnitude of 0(t) for T = 30.0 exceeds that for T = 
at a time t m and t w 2ir. We must notice that the 
lowest value of 0(t) for T = 30.0 located at t = t c « 0.1 
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FIG. 2: The orientation factor 0(t) in quantum paired rotors. 
The figures (a) and (b) correspond to arrangements (a) and 
(b) in Fig. 1. The quantity H/Ek is taken as the unit of 
time. Left: The time dependence for 0(t) within a long time 
scale. The parameter F is varied as V — (dotted line) and 
T = 30.0 (solid line) fixing the kick strength P = 10. Right: 
The behavior of 0(t) within a short time scale. F is varied 
as T = (dotted line), F = 1.0 (dashed-dotted), F = 3.0 
(dashed-dotted-dotted), and F = 30.0 (solid). 



is remarkably smaller than that for T = 0. This indi- 
cates that, in arrangement (a), the orientation of paired 
rotors is enhanced by introducing strong dipolar interac- 
tion. For arrangement (b), in contrast, the magnitude 
of 0{t) for T = 30.0 does not exceed that for Y = at 
any t. In addition, the lowest value of 0(t) at t w 0.1 
seems to be invariant to the change of T, implying that 
the orientation of paired rotors in arrangement (b) is not 
much affected by dipolar interaction. 

In order to examine the effect of the interaction for 
the lowest value of O(t), we investigate in detail the be- 
havior of 0(t) around the focal time t c with varying T. 
The calculated results are shown in the right two images 
in Fig. 2, where the value of V is increased from T = 
(dotted line) up to T = 30.0 (solid line). In case (a), the 
increase in T monotonously reduces the lowest value of 
0(t). For r = 30.0, the factor eventually takes the low- 
est value 0(t c ) = 0.156 at the focal time t c = 9.1 x 10~ 2 , 
which is much smaller than the lowest value of 0(t) for 
r = 0. In case (b), on the other hand, the time depen- 
dence for 0{t) hardly changes with varying T. We thus 
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conclude that, as far as arrangement (a) is concerned, the 
orientation of paired rotors can be efficiently enhanced by 
taking into account the dipolar interaction between ro- 
tors. This is one of main findings of the present study. As 
we see in the next subsection, the enhanced orientation in 
quantum rotors cannot be interpreted from the classical 
dynamics for paired rotors, indicating that the enhanced 
orientation stems from a purely quantum effect. 

B. Classical paired rotors 

Before proceeding to a further investigation of quantum 
paired rotors, we consider the effect of dipolar interaction 
on the orientation for classical paired rotors. For classical 
kicked rotors, the term Hi defined in (2) is rewritten as 

H i = I -6 2 (t) + V(6 i ,t), (16) 

while the interaction term Wij is the same as that de- 
fined in (4). Under field- free conditions, the equations of 
motion for arrangement (a) are given by 

I6i = ; ^[3sin(0 1 +0 2 )-sin(0 1 -0 2 )], (17) 
Ih - ^[3sin(0i+0 2 )+sin(0i-02)]- (18) 

By transforming the variables into q\ = 6\ + 62 and q 2 = 
81 —02, we obtain the following equations: 

q\(t) = 6r c; sin gi (t), (19) 
q 2 (t) = -2T cl smq 2 (t). (20) 

Here we define the parameter T c i = Ed/ (21), showing 
the strength of the dipolar interaction between classical 
rotors. For arrangement (b), the same procedure yields 

q\(t) = 0, (21) 
h(t) = AT cl sinq(t). (22) 

Solutions of Eqs. (19)-(22) are expressed by Jacobi's ellip- 
tic functions. The orientation factor for classical paired 
rotors is calculated by 0(t) = 2 — 2(cosgi(t) cos q 2 (t)) c i, 
where the bracket (• • • } c / means averaging over initial 
angles 9 t (t = 0) [24]. 

Figure 3 shows the time dependence of the orientation 
factor for the classical paired rotors. Quantities I/P and 
(I/P) 2 are taken as the unit of time and the parameter 
r c ;, respectively. The value of T c i is incrementally in- 
creased from T c ; = (dotted line) up to T c i = 45 (solid 
line) as denoted in the figure caption. We see that the 
increase in r c / inevitably raises the minimal value of 0(t) 
in both arrangements (a) and (b). This leads to the con- 
clusion that, in classical systems, the strong dipolar in- 
teraction interferes with the orientation of paired rotors 
for both arrangements (a) and (b). The physical inter- 
pretation is given as follows. When two dipolar rotors 



2.5 1 1 1 . , — r 1 2.5 




0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 
Time f Time f 

(a) (b) 

FIG. 3: The orientation factor 0(t) in classical paired rotors. 
Figures (a) and (b) correspond to arrangements (a) and (b) 
in Fig. 1. The quantity I/P is taken as the unit of time. 
The parameter r c ; = Ed/ (21), determining the strength of 
dipolar interaction, is varied as r c ; = 0, 15, 30, 45 from the 
bottom (solid line) to the top (dashed line). 

are assigned in arrangement (a), strong dipolar interac- 
tion forces them to be parallel in the direction normal 
to the field direction (See Fig. 1). Hence, the interaction 
hinders the orientation of rotors in the field direction. 
For arrangement (b), on the other hand, two dipolar ro- 
tors tend to be anti-parallel to each other. This increases 
the minimal orientation factor at the focal time. As a 
consequence, the dipolar interaction in classical systems 
certainly prevents the rotors from becoming oriented in 
the field direction defined in Fig. 1. 

These facts naturally lead us to the following question: 
Why is it that a strong dipolar interaction can enhance 
the orientation of rotors in quantum systems (See Fig. 2 
(a)) ? Comparing the behavior of 0(t) shown in Fig. 3 (a) 
with that shown in the right of Fig. 2 (a), we clearly see 
that the effect of dipolar interaction on the orientation 
of paired rotors differs completely between classical and 
quantum systems. To settle this point, we consider the 
time development of wavefunctions in the paired-rotor 
system as follows. 



C. Time development of wavefunctions for 
post-kicked time 

To understand the mechanism underlying the en- 
hanced orientation in quantum paired rotors, we ex- 
amine the time development of the probability density 
6 2 , t)\ 2 in the 61-62 space. Figure 4 (i) and (ii) give 
contour plots of the probability density for arrangement 
(a) with T = 30.0. At t = (Fig. 4 (i)), the amplitude of 
\^(6\, # 2 ,t)| 2 is spatially confined around the two sym- 
metric positions (61,62) = (n/2,ir/2) and (— n/2, — ir/2). 
When the kick is applied, these two wave packets move 
toward the origin 6\ = 62 = while retaining their 
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FIG. 4: Contour plots of the probability density \<P(0i,6 2 ,t)\ 2 
for arrangement (a) in Fig. 1. Top: The parameter V = 30.0 
with (i) t = 0, and (ii) t = t c . Bottom: T — 1.0 with (iii) 
t = 0, and (iv) t = t c . White regions indicate large amplitudes 
of the probability density. 



shapes, and finally collide head-on with each other at 
the origin at a focal time t = t c . The resultant angular 
focusing is demonstrated in Fig. 4 (ii), where the prob- 
ability density of \\P(0i, 6*2) | 2 is well localized at around 
the origin. This transient angular focusing at the origin 
enhances the orientation of quantum paired rotors. 

We must note that angular focusing is enhanced only 
when r > 1, namely, when the two rotors are strongly 
correlated via dipolar interaction. If the dipolar interac- 
tion is sufficiently weak (T < 1), the probability density 
for the initial state is broadly distributed in the 61-62 
plane (Fig. 4 (hi)). After a kick is applied, the amplitude 
of the wavefunction spreads out over the 61-62 space, fol- 
lowed by the formation of a "rainbow structure" [12] at a 
focal time t c (Fig- 4 (iv)). The degree of angular focusing 
for the rainbow structure is obviously inferior to that for 
the case of strongly interacting rotors (Fig. 4 (ii)). In 
summary, two factors are essential for enhancing the ori- 
entation of paired rotors: i) The initial state before the 
kick consists of two wave packets strongly confined at 
symmetric positions with respect to the origin, ii) These 
wave packets move translationally toward the origin after 
the kick. It should be mentioned that the translational 
motion of two wave packets after a <5-pulse is not triv- 
ial. We can analytically trace the motion of those wave 
packets by calculating the expansion coefficients Dui ap- 
pearing in Appendix B. Details of the calculations will 



O 




FIG. 5: The orientation factor 0(t) for paired rotors kicked 
with a sequence of seven pulses of the strength P = 10; The 
values of Y are varied as V = 0, 1.0, 3.0, 5.0, 10.0, 30.0 from 
the top (the solid line) to the bottom (the dotted one). The 
quantity h/En is taken as unit of time. 



be published elsewhere [25]. 

The spatial profile of the initial eigenstate is deter- 
mined by the potential term Wu- For arrangement 
(a), the potential W12 ("IS (X function of 6\ and 62 gives 
two potential minima, at (6*1,02) = (7r/2, tt/2) and 
(-7r/2,-7r/2), and a maximum at {6 1 ,6 2 ) = (0,0) [20]. 
The energy difference between the minimum and the 
maximum is determined by the interaction energy Ed 
or, equivalently, the parameter T = Ed/Ek- When V 
is much larger than unity, the energy difference becomes 
so large that the initial eigenstate is strongly localized at 
the two potential minima, as shown in Fig. 4 (i). In addi- 
tion, the amplitude of the wavefunction at the origin be- 
comes almost zero due to the large potential maximum. 
For post-kicked time, however, the kick creates a large 
number of excited states, so that a superposition of them 
can produce a transient angular focusing at the poten- 
tial maximum (61,62) = (0,0). This leads to a minimal 
orientation factor at a focal time. On the other hand, in 
the classical limit, the orientation of paired rotors in the 
field direction (61, 6 2 ) = (0, 0) cannot occur when the en- 
ergy difference between the minimum and the maximum 
is larger than the kinetic energy of the rotors immediately 
after the kick. In other words, the strong dipolar inter- 
action prevents the paired-rotor state from being located 
at the origin (61,62) — (0,0). This follows that the en- 
hanced orientation of paired rotors is a purely quantum 
phenomenon. 



D. The accumulative squeezing 

The orientation of paired rotors can be further en- 
hanced by applying the " accumulative squeezing" scheme 
proposed in Ref.[12]. This strategy is based on a specially 
designed series of short laser pulses leading to a dramatic 
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narrowing of the rotor angular distribution. Figure 5 
shows the orientation factor 0(t) of paired rotors kicked 
by a sequence of seven pulses of the strength P = 10. The 
values of T increase from (solid) to 30.0 (dotted). For 
both arrangement, the strategy works well to achieve the 
angular squeezing in paired rotors. Moreover, in arrange- 
ment (a), a considerable reduction of the factor 0{t) is 
seen for strongly interacting paired rotors with r>l. 
This result provides a new prospect for the scheme of 
multiple-pulse angular squeezing in interacting quantum 
rotors. 



pansion as follows [22]: 

oo 

ce 2n {a,v a ) = ^2 4 2 mK)cos2roa. (Al) 

m=0 

oo 

se 2n+1 {a,v a ) = ^2 B 2l"+i 1) K) sin ( 2m + l)a(A2) 

m=0 

oo 

ce 2 „+i (a,v a ) = £4£+ 1 1) (v Q )cos(2ro+l)a(A3) 

oo 

se 2n+2 {a,v a ) = 4mi 2) K)sin(2m + 2)a.(A4) 



IV. CONCLUSION 



In conclusion, we have theoretically investigated the 
quantum dynamics of paired kicked rotors. The orien- 
tation of paired rotors after the 5-function kick is re- 
markably enhanced by introducing dipolar interaction 
between rotors, when the rotors are deposited in an iden- 
tical plane. The enhanced orientation is attributable 
mainly to two factors: i) The initial state before the kick 
consists of two wave packets strongly confined at symmet- 
ric positions with respect to the origin in the 61-62 space; 
ii) These wave packets move translationally toward the 
origin after the kick. We have also demonstrated that 
the orientation of quantum paired rotors can be further 
enhanced by applying a specially designed sequence of 
pulses. Our findings will stimulate experimental works 
aimed at the orientation of polar molecules correlated 
via dipolar interaction. 



By substituting Eqs. (A1-A4) into Eq. (10), we obtain 
successive relations that determine the expansion coeffi- 
cients. For |^2m'}' as an example, we obtain the fol- 
lowing relation: 



e4 2n) -v4 2n) =0, 



(A5) 

(e - 4)4 2n) - v (2A 2n) - 4 2n) ) = 0, (A6) 

(e 4m 2 )4 2 : i) - v (A^_ 2 A™ a ) = 0. (m (AZ) 

The orthogonality of the Mathieu functions is described 
by 



da cei(v a ,a)cei>(v a ,a) = Tt5 w , (A8) 
da sei(v a ,a)sei>(v a ,a) = nSw , (A9) 

-7T 

/da cei(v a ,a)sei'(v a ,a) = 0. (A10) 
■It 
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APPENDIX A: EXPANSION OF MATHIEU 
FUNCTIONS 



Four kinds of Mathieu functions, ce2 n , se2 n +i, ce2 n +i, 
and se2n+2 , can be expressed in terms of the Fourier ex- 



APPENDIX B: EXPLICIT FORM OF EQ.(14) 

The explicit form of Eq. (14) can be obtained by ex- 
panding the state ^(£,77, t + ) in terms of the Mathieu 
functions. Using these relations, the function (14) is ex- 
panded as 



00 00 



V, r + ) = E D »'f^ V €)91'(V, «,). ( B1 ) 

1=0 Z'=0 

where each type of the Mathieu function is abbreviated 
as 

fv&vt) = ce,(£,^), (B2) 
/2m (Z,v £ ) = se^vz); I = 0,1,2- •• . (B3) 

The definition of gi'(rj, v v ) is the same as that of The 
expansion coefficients {Dui} are calculated straightfor- 
wardly as 

pTT p TT 

0,r = I dU dr, ^,V,r + )fi(Z,v s ) 9ll ( V ,v v ). 



TT J —IT 



(B4) 
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Substituting Eq. (Bl) into Eq. (14), we obtain the ex- 
plicit form of 77, r+ + t) as 



1=0 l'=0 L 

x D H ,fi(£,vt:)gv{r),v v ). (B5) 



In actual calculations for Eq. (B5), the double summation 
with respect to I and I' can be truncated at a finite value, 
because the expansion coefficient Dw rapidly decay with 
increasing I and I'. 
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